Spin-Charge Separation in Two-dimensional Frustrated Quantum Magnets 
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The dynamics of a mobile hole in two-dimensional frustrated quantum magnets is investigated 
by exact diagonalization techniques. Our results provide evidence for spin-charge separation upon 
doping the kagome lattice, a prototype of a spin liquid. In contrast, in the checkerboard lattice, 
a symmetry broken Valence Bond Crystal, a small quasi-particle peak is seen for some crystal 
momenta, a finding interpreted as a restoration of weak holon-spinon confinement. 

PACS numbers: 71.10.Hf,71.20.-b 
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In the last years, there has been growing interest in the 
investigation of frustrated magnets, both on experimen- 
tal and theoretical sides. Spin-1/2 Heisenberg models on 
two and three dimensional (3D) lattices have been stud- 
ied extensively providing increasing evidences for new ex- 
otic phases like the spin liquid Resonating Valence Bond 
(RVB) state or the Valence Bond Crystal {VBC) 13. 
In two dimension (2D), the checkerboard lattice [3j.|_l|.[5[. 
i.e. the 2D analog of the 3D pyrochlore lattice, and 
the kagome lattice 0, (see Fig. , respectively, seem 
to provide simple theoretical realizations of VBC and 
RVB groundstates (GS). In addition, an exotic scenario 
where magnons break up into {S = 1/2) deconfmcd 
spinons 8] could take place in a RVB state. These theo- 
retical conjectures have triggered a renewal of experimen- 
tal activity on quantum frustrated magnets like e.g. the 
(S = 3/2) kagome compound SrCrgpGa^-gpOig |2j and 
the pyrochlore titanates RE2Ti207 (where RE is a rare- 
earth magnetic atom such as gadolinium or terbium) |lfj| . 
These materials exhibit very rich phase diagrams as a 
function of chemical composition, temperature, pressure 
and magnetic field which underline the crucial role of 
geometric frustration. Furthermore, the recent discovery 
of sup er conductivity in a C0O2 based compound with 
triangular layers [ll| has revived interest in exotic RVB 
mechanisms of superconductivity [l^ . 

In this Letter, we investigate the dynamics of a mobile 
hole injected in various 2D frustrated quantum magnets. 
This study is indirectly motivated by the long-standing 
puzzles offered by the "normal phase" of high-T c super- 
conductors whose discovery initiated a search for spin- 
charge separation mechanisms in 2D and the first pro- 
posal of 2D Luttinger liquids 01 • The dynamics of a mo- 
bile hole doped into an (Ising) antiferromagnet (AF) was 
first studied by Brinkman and Rice and a coherent 
motion was later predicted by Trugman However, 
a close inspection of the hole spectral function showed 
evidences for a composite nature of the doped hole quasi- 
particle (QP) showing holon and spinon components (lif . 
Large quantum fluctuations present in frustrated mag- 
nets might then induce a complete deconfinement of the 
hole components. Here we present large scale exact diag- 
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FIG. 1: (a) Schematic picture of the VBC GS of the Heisen- 
berg model on the checkerboard lattice. The two degenerate 
singlet patterns are shown with full blue and dashed red lines, 
(b) Dotted line: Brillouin zone (BZ) of the checkerboard lat- 
tice. Full line: extended BZ equivalent to the square lattice 
BZ. The path from the zone center F to the k = (n, n) point 
M is shown as a dashed line, (c) The kagome lattice; (d) 
Dotted line: BZ of the kagome lattice. Full line: extended BZ 
equivalent to the nondepleted triangular lattice BZ. 



onalization results for single hole spectral functions and 
hole-spin correlations obtained on the checkerboard and 
the kagome antiferromagnets. A finite lifetime of the 
hole (breaking apart into a spin 1/2 spinon and a charge 
Q = e holon) is signaled by the absence of a QP peak in 
the computed spectral functions. Since the latter could, 
in principle, be measured by Angular Resolved Photoe- 
mission Spectroscopy (ARPES) on high-quality (cleaved) 
insulating single crystals |17| , our results can provide mo- 
tivation and theoretical background to experimentalists. 
In addition, this investigation might also be relevant to 
lightly doped samples. Evidences for spin-charge sepa- 
ration in the kagome spin liquid GS are provided while, 
in the checkerboard lattice, we find small QP peaks for 
some momenta. 
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Before going further, let us first summarize briefly the 
properties of the undoped insulating systems we are deal- 
ing with. The checkerboard lattice is believed to form a 
VBC with plaquette singlets on a subset of the void pla- 
quettes 0, 0, 0, ^| . The groundstate manifold is two 
fold degenerate [see Fig. Q] (a)] . The system has quite 
a sizable triplet gap At ~ 0.7 J El lal, and a some- 
what smaller singlet gap As w 0.3 J [j, |ljj to domain 
wall like singlet excitations 19]. Based on this picture 
of a fully gapped system, we expect a coherent QP mo- 
tion of an injected hole at least in parts of the Brillouin 
zone (BZ) as observed for example in Heisenberg ladders 
or dimerized spin chains. On the other hand, the un- 
doped kagome lattice has a very interesting and puzzling 
groundstate of a new type: while this system is also mag- 
netically disordered p, it has quite a small triplet gap 
(Ay ~ 0.05 J) 0- The singlet sector seems to remain 
gapless showing an exponentially large number of sin- 
glets within the spin gap 0, H3 ■ These unconventional 
low lying excitations open the door to new and surprising 
phenomena upon hole doping. 

Our numerical calculations are based on the standard 
t—J model Hamiltonian: 

' E V ( c \,* c w • !;•<••) S < S - 7"<"/ C 1 ) 

{i,j),<r 

where on both lattices all bonds have the same couplings 
t and J. This model is believed to be reliable to de- 
scribe weakly doped Mott- Hubbard insulators with large 
optical gaps. Hereafter, a physical value of J/\t\ = 0.4 
is assumed, unless mentioned otherwise. The spectral 
functions are defined in the standard way: 

(2) 

and calculated by the Lanczos continued-fraction tech- 
nique on finite clusters of up to 32 sites with periodic 
boundary conditions. Here |^o) (of energy E ) is the GS 
of the undoped insulating system, so that the dynamics 
of a single hole is probed. A simple sum rule is satisfied: 
/ A a (k,uj)duj = N a /N ~ 1/2, independent of k. There- 
fore the spectral function plots share an arbitrary, but 
common scale, unless mentioned otherwise. In addition 
an artificial broadening r] — 0.05|i| or 0.1\t\ has been used 
to plot the data. The weights of the lowest peaks can be 
obtained exactly, independent of r\. We have indicated 
this weight in the plots when a QP peak is present. 

As a first step it is of special interest to consider the 
static limit i.e. the case t = (equivalent to an infi- 
nite hole mass) where a spin 1/2 is replaced by an in- 
ert (impurity-like) site. The spectral function, which 
becomes here independent of momentum k, gives useful 
hints on the host magnetic disturbance induced by the 
static hole. Results in Fig. [3 show very different qual- 
itative behaviors for the checkerboard and the kagome 
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FIG. 2: Single hole spectral function for a static (t = 0) hole 
in a checkerboard (a) or a kagome (b) lattice. The checker- 
board spectral function is dominated by the lowest energy 
peak, in contrast to the kagome response. For this case the 
weight of the poles is indicated by red circles. 



lattices. A sharp peak which exhausts most of the spec- 
tral weight is found in the checkerboard lattice. This 
is consistent with the picture of breaking a local pla- 
quette singlet, therefore producing a localized moment, 
which remains confined to the hole. The very large peak 
weight shows that the bare state Cj^l^o) is a very good 
approximation of the true one hole groundstate. This 
however seems to be quite different for the kagome case: 
the spectral function shows a rather broad response, with 
no pronounced low energy peak [see weight distribution 
(circles) in Fig. [2(b)]. Here the overlap of the bare elec- 
tron removal state and the true groundstate of the one 
impurity problem is tiny. This can be nicely understood 
in the "dimer freezing" picture advocated in Ref. [2l| . 
where the screening of the impurity constitutes an im- 
portant deformation of the RVB spin liquid groundstate. 
This suggests that a single dopant in the kagome lattice 
generates a dimer-screened holon with charge Q=e and 
spin 5=0 and a deconfined spinon (Q=0, 5=1/2). 

We now move to the case of a dynamic hole (t ^ 0). 
Because of the absence of particle-hole symmetry in frus- 
trated lattices one has to distinguish between t > and 
t < 0. Note that for t < frustration can also appear in 
the hole motion. For example, a tight-binding particle on 
an isolated triangle gains a kinetic energy |i| per particle, 
a factor of two smaller than for t > 0. Our results for 
a 32 cluster are shown in Fig. [3] along the T <-> M 
line of the BZ [Fig|H(b)]. In all cases most of the spec- 
tral weight is found to be incoherent, distributed over a 
range of 7 — 9 \t\. However, a small QP peak is visible, in 
particular for momenta close to the M -point. The region 
close to the T point has a tiny or no QP peak and the 
shape of the spectral function at T itself is very special, 
probably because of a higher point group symmetry. It 
is interesting to notice, in the case t < 0, the presence of 
a pseudo-gap in the region to ~ 0. 

We have estimated the QP bandwidth W of the 
checkerboard lattice from the dispersion of its QP pole. 
As shown in the left panel of Fig. 21 the single hole band- 
width is much smaller than in a 2-leg spin ladder, which 
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FIG. 3: Single hole spectral functions obtained on a 32 site 
checkerboard cluster along the line F <-> M. Left panel t = 
+1, right panel i = — 1. In both cases J/\t\ — 0.4. When 
a quasiparticle peak is present the corresponding weight is 
indicated. Scaling factors are applied as indicated on the 
plots. 
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FIG. 4: Left panel: QP bandwidth in the checkerboard lat- 
tice vs J/t (J > 0). Comparison is shown with the case of 
the square lattice (data from Ref. flrjl L The magnitude of the 
dispersion of the low-energy edge of the kagome one-hole spec- 
trum is also shown. Right panel: checkerboard QP weight at 
the M point of the BZ vs J/\t\ and square lattice QP weight 
at the E = (ir/2,7r/2) point. 

has a gapped spectrum, (W ~ 2t for strong magnetic 
rung couplings 23]) or even in the square lattice, which 
has a gapless spectrum, where W ~ 2.2J. In the latter 
case, the strong renormalization of the coherent motion 
can be explained e.g. by long-wavelength spin- waves scat- 
tering [jjj UJj ■ In the checkerboard lattice the QP are 
more massive, almost localized, although the VBC host 
has no low energy excitations. We believe that this re- 
markable fact is due to a destructive interference effect 
between the paths available for the hole to hop from one 
plaquette to the next. One can define QP weights Z lh (k) 
as the relative weights contained in the first (5-function 
appearing at the bottom of each spectrum in Fig. [3J Its 
behavior versus J/t is shown in the right panel of Fig. 0] 
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FIG. 5: Single hole spectral functions (black lines) along 
the line F <-> M computed on a 27 site kagome cluster for 
i = +1 (left panel) and t = — 1 (right panel). In both cases 
J/\t\ — 0.4. The red circles denote pole locations and their 
residues. Note that no quasiparticle peaks are visible for all 
momenta. 

and compared with existing data for the square lattice. 
We note that, although there are quantitative differences, 
the qualitative behavior is similar; larger J leads to larger 
holon-spinon confinement and hence to an increase of Z lh . 

The spectral functions of the kagome lattice shown in 
Fig. show exotic behaviors; they are very broad for 
all momenta (width ^6 — 8 |t|) and, in contrast to the 
checkerboard lattice, do not show visible QP peaks both 
for t > (left panel) and t < (right panel). Note that 
the broad appearance is not due to a large ry, but is an 
intrinsic feature of the spectral function as can be seen 
from the large number of poles carrying spectral weight 
(circles in Fig.JSJ). Note that the GS |*P ) of the 27 site un- 
doped kag ome cluster is a spin 1/2 state, e.g. S z = +1/2, 
||). Therefore, the initial state of charge Q = e and spin 
projection S z = 0, Ck.| |T), is an equal- weight super- 
position of a singlet and a triplet state leading to two 
contributions. We have checked that both components 
are consistent with the decay of the original hole into a 
holon and a spinon (with no QP peak). Close similari- 
ties of the singlet and triplet channels also suggest that 
the spinons in the final state (with either parallel or anti- 
parallel spins) are weakly interacting in agreement with 
the mechanism of spinon deconfincment. 

The spectral function data strongly support a spin- 
charge separation scenario for the kagome lattice. In or- 
der to get further evidence and more insight about the 
origin of this phenomenon we also investigate the spin 
density profile in the vicinity of the hole, i.e. the hole- 
spin correlation function (n^Sf) in the one hole-doped 
30 site kagome cluster (for which Sz = 1/2). Results are 
shown in Fig.EJfor t — +1 (a) and t = — 1 (b). In the case 
t = +1 we clearly see a repulsion between the net 5=1/2 
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FIG. 6: (a-b) Hole-spin density-density correlation function 
{n^Sf} calculated on a 30 site kagome cluster for both signs of 
t. The diameter is proportional to the magnitude of the cor- 
relation. Full blue (striped red) circles denote positive (nega- 
tive) correlations, (c-d) Nearest neighbor bond spin-spin cor- 
relations around the mobile hole computed on a 27 site sample 
for both signs of t. J = 0.4 for all the plots. Full blue lines 
denote antiferromagnetic correlations, while the red, dashed 
lines stand for ferromagnetic correlations. 

moment and the mobile hole. The case t = — 1 displays 
a different picture, showing a rather uncorrelated behav- 
ior between the spinon and the holon. Nevertheless a 
spinon-holon confinement does not seem to be present, 
in agreement with spectral function data. In contrast 
we have checked that on a checkerboard lattice a clear 
localization of the spinon close to the holon takes place. 

Finally we consider the GS of the one hole-doped 27 
site kagome cluster of total charge Q = e and spin 5 = 
to mimic the holon wavefunction. Figs. (c,d) show a 
snapshot of the nearest neighbor spin correlations around 
the mobile holon. Surprisingly in the case t = +1 a 
"dimcr freezing" is observed very similar to the static 
impurities considered in Ref. |2l|. In contrast, the case 
t = — 1 shows two slightly ferromagnetic bonds on the 
two triangles containing the holon, while the other bonds 
remain largely unaffected. This behavior of the magnetic 
environment is fairly well described by a simple energetic 
consideration on a 3-site ring t—J model; for t = +1 the 
propagating holon is accompanied by a singlet, while in 
the opposite t — —1 case an accompanying ferromagnetic 
bond appears. 

To conclude, based on extensive analysis of the single 
particle spectral functions and various charge-spin corre- 
lations in the single hole GS, we have addressed the issue 
of spin-charge separation in 2D magnetically disordered 
magnets. We show strong evidence that it occurs in the 
kagome lattice and we expect it to be robust for small but 



finite doping. This provides the first example of observed 
spin-charge separation in a two-dimensional microscopic 
model. The RVB nature of the undoped groundstate 
(characterized by a gapless low-excitation singlet spec- 
trum) seems to be crucial for this behavior. Indeed, in 
the checkerboard lattice, which exhibits a VBC struc- 
ture, a weak hole-spinon confinement manifests itself as 
QP peaks for some momenta. 
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